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Abstract We investigate phase-sensitive interference effects in a periodically sin(2-7r/ r fi)- 
driven, artificial two-state system connected to a microwave resonator at /r,c — 800 
MHz. We observe two kinds of multiphoton transitions in the two-state system, accom- 
panied by: 1) Several quanta from the drive at / r f and 2) one quantum at / r f and several 
at /lc- The former are described using phase-sensitive Landau-Zener transitions, while 
the latter are discussed in terms of vibronic transitions in diatomic molecules. Inter- 
ference effects in the vibronic transitions governed by Franck- Condon coefficients are 
also considered. 

Keywords Landau-Zener tunneling ■ superconducting qubits ■ multiphoton 
transitions ■ Franck-Condon physics ■ Aharonov-Anandan phase 



1 Introduction 

Interference effects have been found to play an important role in the dynamics of qubits 
1,2,3,4,5,6 . For their understanding, the interference phenomena can be considered 
from different points of view: as interference between successive Landau-Zener (LZ) 
tunneling events [7][8][9], analogous to Mach-Zehnder interference [10], or in terms of 
spin dynamics |11] . In charge-phase qubits, strongly coupled to a microwave resonator, 
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analogy with vibrational transitions in diatomic molecules has been pointed out [12] . 
The fact that the latter phenomena are related to interference effects has its origin in 
the phase-space dynamics |13| . 

The state of a superconducting charge-phase qubit or a Cooper-pair box (CPB) 
can be monitored continuously by measuring its reactive response, either as an effective 
capacitance or inductance, which both are due to the curvature of the qubit's energy 
bands with respect to charge or phase degrees of freedom. When such a system is 
made as a part of an electric LC-oscillator circuit, any change of reactance, caused 
by the evolution of occupancies of the qubit's energy levels, will shift the resonance 
frequency of the resonator. This scheme has been used in several recent experiments, 
typically having the qubit, with splitting AE, connected to a coplanar stripline cavity 
near resonance |14] . We have, however, worked in the fully detuned limit where fag <C 
AE/h. At weak drive, such a read-out scheme perturbs the investigated system only 
weakly, allowing studies of inherent interference phenomena of the periodically driven 
two-level system. 

We have investigated interference effects in the Cooper-pair box and in a charge- 
phase qubit circuit, configured as dual to the CPB (i.e., connected to the resonator via 
the phase, rather than the charge port), see Fig. [1] Both circuits provide good model 
systems of interference effects in periodically driven qubits. In the latter case, a clearly 
stronger coupling between the qubit and the resonator could be achieved, which made 
a difference in the observed effects. 

The interference effects in these two circuits are diverse mostly due to the differ- 
ence in the coupling between the resonator and the qubit. On the whole, interference 
phenomena can be understood in terms of multiphoton transitions. In strongly driven 
qubits, clear multiphoton transitions have been observed when the energy quantum of 
the rf-drive, multiplied by a small integer, matches the qubit level splitting 15 , 16 . In 
fact, multiphoton transitions allow a clear-cut distinction between our main findings. 
We may classify our observations as transitions where: 1) Several quanta from the drive 
at / r f and 2) one quantum at / r f and several at f^Q are exchanged between the qubit 
and its surrounding entities, i.e. the classical rf-drive port and the LC resonator. The 
first class can be treated as a sequence of phase-sensitive Landau-Zener transitions 
which produce a characteristic interference pattern, dependent on the "finesse" of the 
interferometer, namely the decoherence of the system. The latter class of transitions 
can be described in terms of vibronic transitions in diatomic molecules where potential 
changes are fast compared to the vibrational frequencies and non-adiabatic coupling 
between levels differing by a large number of quanta becomes possible. These phenom- 
ena include creation of several resonator quanta from the ground state, which is the 
non-adiabatic behavior referred to as the dynamical Casimir effect [17] . 

In this paper, we summarize some of our recent results 4,12 on interference ef- 
fects in a driven superconducting charge qubit coupled to a resonator. The original 
references, in particular, used theoretical considerations and numerical simulations of 
the Bloch equations and the linear-response approximation for a detailed analysis of 
the shape of interference fringes when dissipation is important. Instead, here we try to 
concentrate on robust features which are amenable to general conclusions. For example, 
part of our analysis will be dealing with the phenomenon of destruction of tunneling. 

The paper is organized as follows: First, we start with basics of the Cooper-pair 
box in Sect. 2 and describe the dynamics of a two-level system in Sect. [3] paying 
special attention to the destruction of tunneling from the ground state to the excited 
state. Sect. [4] covers our experimental work and results on CPBs, and it includes some 
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Fig. 1 Schematics of our experimental configurations: a) Cooper-pair box and b) the induc- 
tively read charge-phase qubit (LSET). In both configurations the qubit-induced change in the 
"quantum" reactance alters the resonant frequency £j£c/27r of the tank circuit. This change 
causes a shift in phase of the reflected, constant-frequency microwave signal as illustrated in 
c). In the illustration, C e g grows which leads to a reduction of w^c and, consequently, to a 
decrease of the reflection phase arg(i~'). 



considerations aiming towards understanding of the Stokes phase and its relation to 
geometric phases in a spin- 1/2 system. We present data at phase bias n (minimum 
band gap), and we analyze, especially, data in the regime of destruction of tunneling 
at the charge degeneracy. Results at zero phase bias can be found in Ref. [4]. Section 5 
describes the analogy of spin flips in strongly coupled qubits with vibronic transitions 
in diatomic molecules. We do not describe much our experimental techniques, but to 
make the paper self-contained, we have explained some of the main technical points 
when presenting the experimental results. 



2 Cooper-pair box 



A split Cooper-pair box (CPB) is formed by a single-Cooper-pair transistor (SCPT) 
embedded into a small superconducting loop 18 ,19 20 , see Fig. [1] The charging energy 
of the CPB, E c = e 2 /2C s ~ 1 Kelvin, is given by the total capacitance C^j which 
includes the junction capacitances, the gate capacitance C g , and the self-capacitance 
of the island. The effective Josephson energy is given by the sum of the energies of the 
individual junctions (-Eji + -Ej2) cos(0/2) = Ej cos(0/2), which is tunable by magnetic 
flux & , i.e. by the superconducting phase across the two junctions, <j> — 2n$/<Po- Here, 
^0 = h/2e is the superconducting magnetic flux quantum. 
We may write the Hamiltonian in the charge basis as 



E c (n - n s ) 2 \n)(n\ - ^ cos ± (\n - 2>(n| + |r> + 2}<n|) 



+^sin^(|n-2)<n|-|n + 2)(ni) 



(1) 
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Fig. 2 The two lowest energy bands of Cooper-pair box as a function of the gate charge n g 
(in units of e). In sweeping n g (horizontal arrows), Landau-Zener tunneling between the bands 
takes place close to the point where the energy difference is at minimum (curved arrows). 
Successive sweeps can lead to interference of the tunneling events. The interference depends 
on the phase JBJ accumulated in the shaded area. 



Here n denotes the number of extra electron charges on the island, and n g = CgVg/e 
is the charge in electron units induced by the gate voltage Vg on the gate capacitor 
with capacitance C g . ii is conjugate to 9/2, where 9 is the superconducting phase on 
the island. The asymmetry of the two Josephson junctions of the CPB is described by 
d= (E n -E 32 )/E 3 . 

By assuming that E c 2> Ej and that n g is close to 1, one can reduce the circuit to 
a two-state system, with the Hamiltonian 

a If -4LE c (l-n g ) -£j [cos (0/2) -Md sin (0/2)] 

2 \ -E 3 [cos (0/2) - id sin (0/2)] 4E c (l-n g ) 

in the basis of the relevant charge states |0) and |2). When this energy operator is 
expressed as a linear combination of the Pauli matrices// = — \(BzOz + B x <r x +ByOy), 
we get the magnetic field components as B z = 4Z? C (1 — n g ), B x = -Ej cos (0/2) and 
By = dEj sin (<f>/2). The energies of the two states as a function of n g are illustrated 
in Fig. [3 

In the experiments Ej — 0.6 K and d = 0.22, which means that the magnitude for 
the off-diagonal components \A\ ranged over 1.4-6.2 GHz [4]. The non-linear paramet- 
ric capacitance, which has been the cornerstone of our interference studies, has been 
experimentally studied in Refs. 521,22 . 



3 Landau-Zener interference 

We study the effect of sweeping the gate charge n g . We denote the energy eigenstates 
by |— ) and |+), and suppose the system is initially in the lower energy state |— ), see 
Fig. [2] Two cases can be distinguished. When the sweep rate is small compared to 
the energy level difference, the system stays in the lower state. The alternative case 
is that the sweep rate is comparable to the level spacing. This may take place close 
to the degeneracy points, where the energy bands would cross in the absence of the 
Josephson coupling. There the system can tunnel from the lower state to the upper 
in a process known as Landau-Zener tunneling. The tunneling probability in a single 
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sweep is given by [7ll8ll9llTT 



> -2tt7 2n A 



Here v is the speed at which the sweep passes the crossing point v = |d(eo — 62)/dt 
and eo — e 2 = 4_B C (1 — n g ) is the energy difference in the absence of the Josephson 
coupling. In Eq. (|3} the adiabaticity parameter 7 defines whether the development is 
adiabatic (7 ^> I) or sudden (7 <C I). 

When the system is coherent, and the degeneracy point is crossed several times, 
the transition amplitudes for each subsequent pass have to be added for evaluating the 
transition probability. These amplitudes may be tracked by the 'scattering' matrix S, 
defined by (cf. Refs. 



HA /Vl-^Lzexp(^ s ) 

\+)J V VT 



(4) 



where phase factors are chosen to simplify the matrix. Here the scattering phase is 
05 — 05 — 7r /2, where the Stokes phase <j>g depends on the adiabaticity parameter 7: 

= tt/4 + arg(r(l - 17)) + 7(ln 7 - 1). (5) 

In the adiabatic limit, 4>s — > — 7r/2, while in the sudden limit, ^5 = — 7r/4. 

Away from the degeneracy point, the eigenstates | — ) and +) accumulate the rela- 
tive dynamical phase 

Vd = - ¥> (_) =~f [E+{n E {t)) - £-(n g (i))] « • (6) 

Thus, the condition for constructive interference is that i^l — an d VR — 2<^g are 
multiples of 27r, where ipi, and (^r refer to the dynamical phases accumulated on the 
left and right sides of the degeneracy point, respectively. For example, in the adiabatic 
limit, i£l,R have to be odd multiples of tv. 

The LZ interference can also be employed to suppress the tunneling to the upper 
level 26 . In fact, the contrast for destructive interference in the experimental data looks 
often better than for constructive interference, and these conditions can be employed 
more simply to determine the behavior of the system as the ground-state response in 
the measurement is well known. In this case the interference conditions are 

ip L + 2nl = 2(f) s and ip R + 2n£ = 2(p s (7) 

with an integer I. 

In the sudden limit, 7 <C 1, the periodically driven two-level system can be solved 
in a rather straightforward manner |27j . At degeneracy, the probability of being in the 
upper state varies periodically with time according to 

P(t) = sin 2 [B x J {A/riu T f)At/h] , (8) 

where Jo denotes the zeroth Bessel function, and A (= 4E c rig below) is the drive 
amplitude at frequency / r f = w r f/27r. One notices that the destruction of tunneling 
takes place, when the ratio A/huirf coincides with a zero of Jo- 

The description of interference patterns in a Cooper-pair box using the scattering 
matrix formalism has been discussed in Refs. [281129] . In general, relaxation phenomena 
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should be included as they have strong influence on the sharpness of the interference 
fringes. For this reason, our main analysis method has been based on spin-1/2 NMR 
simulations [4ll29| . which also includes linear response calculations in order to obtain 
the measured, effective capacitance. Alternatively, the dressed-state approach may be 
employed as has been done by C. Wilson el al. [6]. 



4 Measurement results on CPB 

We have performed low-dissipation microwave reflection measurements 30 ,31,32 on a 
series LC resonator in which the box effective capacitance, 

± _ 8 2 E±(<p,n g ) _ Cl d 2 E±{<j > ,n s ) 
cff 8Vi e 2 dn\ 

is a part of the total capacitance C§ + C where the superscript ± refers to ground 
and exited states of the qubitQ The resonator is formed by a surface mount inductor 
of L = 160 nH. With a stray capacitance of Cg = 250 fF due to the fairly big lumped 
resonator, the resonant frequency is /o = 800 MHz, and the quality factor Q ~ 16 
is limited by the external Zq =50 Q.. When varies, the phase and amplitude of 
the reflected signal V ou t = rV\ n change, which is measured by the reflection coefficient 
r — (Z — Zq) /(Z + Zq) = -Toe 1 arg ( r ) . Here, Z is the resonator impedance as seen from 
the end of the 50 fi coaxial cable used for the reflection measurement. The variation 
in arg(r) due to modulation in C~L is up to 40° in our measurements, corresponding 
to a shift of resonance frequency Af p ~ 6 MHz. In all the measurements, the weak 
probing signal V- ln at frequency f m was continuously applied, in addition to the DC- 
bias and the rf-drive. Thus, the total gate charge variation, in units of e, can be written 
as ng(t) = n g o + n g f sm(2n f^t) + fig 1 sin(27r/ m £), where the I s *, 2 nd , and 3 rd terms 
correspond to the DC, rf, and measurement drives, respectively. 

We have made extensive scans of the CPB reflection by varying the LZ drive fre- 
quency / r f = 0.1-20 GHz and its amplitude n g f — 0-3 electrons, as well as the qubit 
DC-bias: n g o and <j>. In Fig. [3] we present the reflection phase arg(_T) and magnitude 
\r\ measured at / r f = 4 GHz on the plane spanned by n g f and n g Q. The Joseph- 
son capacitance C c g ~ — arg(_T)C^ 2 Zo/(2\/i) deduced from the data of Fig. [3] has 
been given in Ref. ^Jf|. We observe in Fig. [3] a clear interference pattern whose main 
features confirm the coherent LZ tunneling picture: 1) There is an onset of the interfer- 
ence speckles, where the rf drive just reaches the avoided crossing (charge degeneracy), 
n g o ± n g f = n gi dcg = °dd integer, with a linear dependence between n g o and the AC 
drive amplitude. Additionally we have shown in Ref. |3] that 2) the density of the dots 
is proportional to l// r f in the direction of n g o as well as n g , and 3) the pattern loses 
its contrast below a certain drive frequency, in this circuit around / r f ~ 2 GHz. 

In Fig. [5^, at charge degeneracy, we observe a clear sequence of blue dots that signi- 
fies a similar response as in the ground state. These dots correspond to the destruction 
of tunneling due to periodic rf-drive of our two level system. Note also that these lo- 
cations correspond to a large reflection magnitude in Fig. [3Jd. In fact, the reflection is 
enhanced compared with the undriven level, which indicates transfer of energy from 

1 For a more detailed description of C c ff, see Refs. 21 ,4 

2 Full circuit analysis was employed in the evaluation of the reported effective capacitance. 
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Fig. 3 Reflection phase (left-a) and magnitude (right-b) measured in the configuration of Fig. 
la at rf-drive frequency of f T { = 4 GHz using phase bias <j> = it, corresponding to the level 
repulsion of 2A = Ejd = 2.7 GHz. n T J specifies the amplitude of the rf drive in electrons and 
n g o is the DC charge bias. The vertical bars give the color scales in degrees and dB for the left 
and right frames, respectively. 



the microwave drive, via the qubit, to the resonator. Apparently, the destruction of 
tunneling prevents the microwave energy from being deposited to the two-level system, 
and the 'extra' energy is dumped out at both drive frequencies. 

The order number of the destructive interference dots in Fig. [3] is plotted in Fig. 
21 At large rf-drives, the interference dots display a linear dependence on the drive 
amplitude. The linear dependence is a sign of asymptotic behavior that can be obtained 
either from Eq. ((7} or from Eq. (|8]). Starting from Eq. (O, we may evaluate </>l,R under 
a drive of AE = 4:E c n g sin(w r fi). This yields 



8E c nl { 
<PL,R = r - 



where m indexes the fringe number, and — i-^ — m, where <?l R are the integers from 
Eq. ©. Thus, we find 



(9) 



(note that the order of resonance grows with <?l, an d Eq. © holds, of course, only 
for sufficiently large ^Li when it gives rig f < 0). For the destructive interference at 
charge degeneracy in the sudden limit this yields rig = (2irhu> T {/8E c ) (<*l ~ j) ■ By 
fitting Eq. ([9} to three/four highest drive points of the data in Fig. [4] we get tps/ 11 — 
0.30±0.03, which is close to the theoretical expectation of 0.25. However, there is always 
an additional contribution from the dynamical phase picked up near the degeneracy 
point where the bands are curving from linear. This contribution was estimated to be 
about 1/3 from the expected Stokes phase. This correction will drop the vertical offset, 
and reduce the Stokes phase down to 0.207T. Hence, taking all the uncertainties into 
account, our final estimate for the Stokes phase becomes (0.2 ± 0.05)-7r in the sudden 
limit. 

A fit using the plain dynamic phase of Eq. ((6| , and the two level approximation for 
the adiabatic energy levels, is given by the dashed curve in Fig. [4] The z-scale of the 
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(arb. units) 

Fig. 4 Order number of destructive interference fringes at charge degeneracy (n g o = 1) as 
a function of the microwave drive amplitude n 1 ^ for phase bias = 0. For comparison, the 
dynamical phase calculated from Eq. (0 is given by the solid and dashed curves with and 
without Stokes phase, respectively. The zeroes of Bessel function Jo are given by the open 
circles. For details, see text. 



calculated curve has been adjusted to match the experimental results at large drive 
amplitudes. The solid curve, on the other hand, is the dynamical phase corrected using 
the calculated Stokes phase in the sudden limit. The curve agrees with data except 
at the lowest points where a possible inaccuracy in the asymmetry parameter has its 
strongest influence. The open circles display the result of Eq. {§), i.e. the zeroes of Jo, 
which are also seen to coincide well with the data. 

The Stokes phase is related to the non-adiabatic geometric phase, the Aharonov- 
Anandan phase 33 ,24 . It is rather straightforward to show that 

*gcom = 2(1 - P LZ )<PS + tt(1 -m) + 7vP LZ (2m + 2£ L - 1). (10) 

From this relation one gets at the charge degeneracy (m = 0): <Pg C om = (l — PLz)(2 i PS~~ 
7t(2^l — 1)), which is nearly the same relation as what we derived for the asymptotic 
fit in Fig. |4j cf. Eq. ((9} • Consequently, the geometric Aharonov-Anandan phase could 
be approximately determined by taking the measured drive amplitude at the point 
of destructive interference and multiplying the result by — (8E c /hu} r f)(l — Plz)- For 
the data in Fig. [4] at 4 GHz, we get a nearly constant value (0.20-0.21) x2-7r at drive 
amplitudes around n g = 0.5-1. 



5 Artificial molecule 

In the inductive SET, LSET, of Fig. QJ>, a stronger coupling between the qubit and 
the resonator can be achieved than in the Cooper-pair-box configuration. This leads 
to more involved multiphonon phenomena which are reminiscent of the transitions in 
diatomic molecules. These transitions can also be viewed as interference effects, and it 
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is convenient and instructive to describe them as interference in the phase space [131 

Ell- 
in order to understand the molecular analogy, we consider the single-Cooper-pair 
transistor as a two-state system. The two states are the analog of two electronic states 
in a diatomic molecule. The transistor is coupled in parallel with an LC oscillator. The 
LC oscillator is the analog of nuclear vibrations in the molecule. Because of the coupling 
between the two-state system and the oscillator, a transition between the electronic 
levels is often accompanied by a change of the vibrational state. Such transitions, where 
both vibrational and electronic quantum numbers change simultaneously, are known 
as vibronic. 

For quantitative analysis, we write the Hamiltonian in the basis of two relevant 
charge states of the qubit: 

H (*, q) = \ (a,E A - a x Ej cos g + a^Ejdsin + ^ + ^ ■ (11) 

Here the flux # and the charge q in the LC oscillator are conjugate variables, and 
(p x ,Oy,o z) are the Pauli matrices. The capacitive energy E c \ = 2e 2 (n g o — 1)/C"x; can 
be controlled by the gate voltage. Another control parameter is the flux bias $b through 
the loop containing SCPT and the inductor. 

In the following we concentrate on the limit, where the oscillator frequency 1/s/LC 
is much lower than the qubit level difference. This means that all changes in the qubit 
system are much faster than in the oscillator system. Therefore, one can diagonalize 
the qubit part separately from the oscillator. The two energies of the qubit are obtained 
by diagonalizing the 2x2 matrix part of Eq. pip : 

E = ±-Je 2 , + E 2 cos 2 P + E 2 d 2 sin 2 P . (12) 
2 V cl J $o J $0 

This form is valid when the Josephson coupling is a small perturbation compared to 
the capacitive energy _E c i. It is also possible to consider the general case of an arbitrary 
ratio of the Josephson coupling to E c \, in which case the energies of the two lowest 
states are given by Mathieu characteristics. This means that the dependence of the 
eigenenergies on $ would be more complicated than in Eq. (|12p but otherwise the 
following analysis remains intact. 

Considering now the oscillator, we can think about the capacitive and inductive 
terms in the Hamiltonian as the kinetic and potential energies. In addition to the 
inductive potential, there is a potential arising from the qubit, since the qubit energy 
(|12l) depends on the flux <5. Thus there are two potential curves for the oscillations 

v±m . + El ^ + £j v * - + ttzi£. , 13 , 

Such a potential is illustrated in Fig. [5K- The different slopes of the qubit energies give 
rise to a relative shift of the minima of the oscillator potentials. The difference in the 
curvatures of the qubit energies shifts the vibrational frequencies. Figure [5^,, which is 
drawn to scale with realistic parameters of our circuit, shows that vibronic transitions 
can be induced by microwave radiation. 

The intensities of vibronic transitions are traditionally analyzed in terms of the 
Franck-Condon principle, and this analysis turns out to be useful for understanding of 
our data. In its classical form, the Franck-Condon principle says that transitions are 
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Fig. 5 The potentials and phase-space picture of vibronic transitions, (a) The potentials U± 
(solid curves) for vibrations correspond to two qubit states ±. The qubit energy il 1 2[) is given 
by dashed lines. The energies of six lowest vibrational states corresponding to both qubit states 
are drawn. The bars represent energy shifts induced by high and low frequency drives, (b) The 
overlap of +,n) states with n = 0, 1, . . . ,7 vibrational quanta with state | — ,3) in the phase 
space. The dots show the centers of the ellipses corresponding to qubit states ±, which are 
displaced from each other by A<P. 

possible between vibrational states, the trajectories of which intersect in phase space, 
and they are most intense between states with coincident turning points. This is based 
on the idea, already stated, that the transitions between electronic states are fast and 
the vibrations are slow, so that the vibrational coordinates, here the flux and the 
charge q, have no time to change during the transition. In addition, &(t) is slowest near 
the turning points, and therefore transitions there are most likely. In the following we 
discuss how this picture can be extended by semiclassical analysis in phase space, as 
shown in Refs. [13134] , and how this interpretation is related to the data obtained. 

The phase space formed by the coordinate and the canonical momentum q is 
illustrated in Fig.[5jD. The vibrational states correspond to elliptic rings. By scaling the 
coordinates properly, a state with n vibrational quanta in a harmonic oscillator can be 
associated with a circular ring between radii y/n and \/n + 1 (Planck-Bohr-Sommerfeld 
bands [34]). This is the case for the upper qubit state in Fig. [5p. Due to a different 
resonance frequency, the vibrational states corresponding to the lower qubit state then 
appear as slightly squeezed in the q direction. The shift of the minima of the potentials 
U± (Fig. (5ji) appears in the phase space so that the ellipses corresponding to the states 
\ + ,n), where + means the upper qubit state, have different centers than the ellipses 
corresponding to the states | — ,n). Vibronic transitions can occur only between states 
whose ellipses overlap in the phase space. In the example of Fig. [5}d, the state — , 3) 
has overlap only with states with 1 < n < 6. In general, the transition rates are 

determined by the area of the intersection of the respective bands. This area is largest 
when the ellipses touch each other tangentially, i.e., when the turning points coincide, 
and this leads to the classical version of the Franck-Condon principle. In the present 
case though, the overlap consists of two crossings of the ellipses, and one has to add 
the transition amplitudes taking into account the relative phase accumulated between 
the crossings. The matrix element of the transitions, (+, n\ — , m) | 2 , is given by the 
formula [13U34] 



(14) 
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Fig. 6 Measured amplitude \T\ of the 
reflection coefficient in the bias plane 
of gate charge n g o and phase (p = 
2-7r^ b /^o- Several concentric circular 
fringes are visible below the V-shaped 
light region. Out of them only the pure 
electronic transition (n = m) remains 
visible at <j> = 7r, and corresponds qubit 
level separation at / r f = AE/h = 22 
GHz. The color bar gives the scale for 
the magnitude of the reflection coeffi- 
cient in dB. The measurement power 
was Pa = —129 dBm referred to the 
coupling capacitor. 
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Fig. 7 Variation of reflection coeffi- 
cient along three fringes in Fig. [6] They 
correspond to n — m = 0, -1, and -2 and 
are denoted by o, •, and A, respectively. 
The fringe behavior is fitted to Bessel 
functions Jq, jf, and j|. For details, 
see text. 



Here Snm is the area of one of the crossing of the elliptic rings, which sets the maximum 
probability of the transitions. The phase in the phase factor is determined by the 
area A nm between the two alternative paths between the crossing points. Eq. (|14|l 
reproduces well the results of full quantum calculation except near the classical turning 
points (where the two crossing areas merge). In the case of Fig. [S] the oscillations are 
nearly harmonic and the relative difference of frequencies is small. When, in addition, 
the difference in the vibrational quanta is small, |n — m\ <C n, the probability P n m 
scales as Jn-mi^by/n) , where b = \/ujC ' /2hA(S> is the dimensionless separation of the 
potential minima. Owing to the parabolic bands in the region near the minimum gap, 
b oc ($ h /<P - 1/2), when \<P h /$ ~ V 2 | < 1- 

The measurements on the artificial molecule are made by analyzing the reflection 
of microwaves at a frequency fo close to the resonance frequency f^c of the circuit. In 
order to induce vibronic transitions, a microwave excitation at frequency / r f was used. 
The resonance condition for vibronic transition between states |— ,m) and |+, n) is 

/rf ^ — + (n-m)f LC , (15) 

where AE is the qubit energy splitting. A measurement of reflection coefficient T in 
the bias plane (n g ,4>^) is shown in Fig. [6] The vibronic transitions are seen as fringes, 
making half-circles in the bias plane. The fringes are located within a crescent-shaped 
area. The language of interference in phase space allows one to easily understand, for 
instance, the shape of this area, the positions of the fringes and the modulation along 
the fringes. At ^ = <?o/2 only the pure electronic transition is seen. This agrees with 
the analysis above since there the displacement of the oscillation minima vanishes. With 
increasing $ ^ more fringes appear, which is in agreement with the growing separation of 
the minima, and the range of the fringes at a given 4>^ can be determined as in Fig.[5]3. 
In addition, the above picture of enhanced probability of the extremal trajectories is 
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reflected in the shape of the absorption pattern in Fig. [6] where the edges of the fan- like 
structure are rather clearly expressed. 

The strength of the fringes in the reflection measurement depends on two separate 
factors. One is the transition probability to the upper qubit state. In large part of the 
bias plane (below the the V-shaped light region in Fig. [6jl the averaged frequency of 
being part time in the upper and part time in the lower qubit state matches better to 
the measuring frequency than being in the lower state only. This leads to absorption 
which is approximately proportional to the transition rate to the upper qubit state. 
The other factor is the heating or cooling effect caused by the high-frequency / r f 
radiation on the vibrations. For a vibronic transition with n > m (i.e. / r f > AE/h), 
the resonance deposits energy into the resonator causing heating of the vibrations. This 
leads to decreased absorption of the measurement wave. In the opposite case n < m, 
the absorption of a high frequency photon cools the oscillator by reducing its quanta 
which leads to increased absorption of the measurement wave. This asymmetry of the 
fringes is clearly visible in Fig. [6] where the cooling fringes at larger radii are stronger 
than the heating fringes at smaller radii. 

Numerical simulations of the artificial molecule has been done in Ref. [12] ■ These 
are based on Bloch equations describing the qubit and classical equations describing the 
circuit. They reproduce well the observed vibronic spectrum, both the matrix elements 
and the heating/cooling asymmetry. Here we present a simplified analysis by comparing 
three of the fringes with the transition probability P nm = J„_ m (26\/n), neglecting the 
heating/cooling effect. Fig. [7] displays the measured reflection magnitude along fringes 
with n — m = 0, -1, and -2 as a function of the flux bias of the qubit. Squared Bessel 
functions Jq, Jf , and jf have been fitted to the data. For the a;-scale argument we 
have taken const x (^b/^0 — 1/2) although this approximation will be rather crude 
when > 0.54. Nevertheless, the comparison can be employed to look for the 

presence of basic interference phenomena. The agreement between the data and the 
fits is quite good, though the modulation of the pure electronic fringe is a bit weaker 
in the measurement than given by Jq. 

The reason for smearing of the Bessel modulation is the variation of the number 
of quanta in the resonator. Assuming the resonator is in a coherent state, it is a 
superposition of the number states and this should be taken into account in the analysis 
of the interference fringes. This diminishes the sharpest features of the interference, but 
does not remove it. Consequently, we may conclude that the dip at ~ 0.54 in the 

fringe magnitude corresponding to the pure electronic transition, is due to destructive 
interference in phase space in formula Eq. (|14[l. 

In summary, we have discussed interband transitions in superconducting Cooper- 
pair boxes and charge-phase qubits and argued, how interference effects in various 
forms can be found to underlie the observed phenomena. 
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